Accurate ionization potential of gold anionic clusters from density functional theory 

and many-body perturbation theory 
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We present a theoretical study of the ionization potential in small anionic gold clusters, using 
density functional theory, with and without exact-exchange, and many body perturbation theory, 
namely the Go Wo approach. We find that Go Wo is the best approach and correctly describes the 
first ionization potential with an accuracy of about 0.1 eV. 

PACS numbers: 31.15.xm,31.15.E-,,33.15.Ry 



INTRODUCTION 

The peculiar properties exhibited by gold clusters 
made them a popular subject for physical chemistry stud- 
ies and numerous applications in recent years, [l|~l8|. Un- 
like bulk gold, which is essentially chemically inert |9|, 
gold clusters possess edge atoms that have a lower coor- 
dination and thus a higher reactivity. As a consequence, 
a close interplay between the structural and electronic 
properties of different Au nano-particles is observed and 
the determination of clusters sizes and shapes has a fun- 
damental importance for both basic studies and applica- 
tions. 

Presently, the atomic structure of gold clusters can be 
established by a number of different experimental tech- 
niques [10| including photoelectron spectroscopy (PES) 
LMI, ic '""' 



; phott 
[l3,i 



nHlSf, ion mobility [iJI, infrared spectroscopy [20|,l2l|, 



trapped ion electron diffraction [22l - l25j . and photodis- 
sociation 26l428|. However, in all cases detailed infor- 
mation from theoretical investigations must complement 
the experimental data to allow the correct interpreta- 
tion of the results and the final resolution of the different 
isomers. In this context, PES is a typical case, where 
photoelectron spectra obtained from theoretical calcula- 
tions on a set of plausible structures are compared to the 
experimental peaks to assign the different features and 
identify the structure of the clusters. 

Usually, the photoelectron spectrum of gold clusters 
is simulated by density functional theory (DFT) calcu- 
lations, since these are easily affordable for a large va- 
riety of structures and provide direct access to approx- 
imate photoelectron spectra through the inspection of 
single-particle orbital energies (within the Kohn-Sham 
(KS) scheme [29|). However, generally the DFT calcu- 
lations are performed using the generalized gradient ap- 
proximation (GGA) for the exchange-correlation (XC) 
functional, and thus several drawbacks are introduced. 
The most important in the present context is the limita- 
tion regarding the description of the single-particle spec- 



trum by these functionals. In fact, correct values of the 
orbital energies cannot be achieved in DFT when an ap- 
proximate description of the exchange interaction is em- 
ployed j30l432|. Due to the self-interaction error (SIE) 
33l |34| present in GGA functionals, orbital energies are 
strongly overestimated and for a comparison with exper- 
imental data a rescaling of the energies is required based 
on additional calculations and/or experimental evidence 
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. Alternatively, for an accurate estimation of the 
first ionization potential (IP), ASCF calculations, based 
on total energy differences between the anionic and neu- 
tral species, can be considered [36|. 

More accurate orbital energies can be achieved with 
exact exchange (EXX) KS approaches, i.e. the opti- 
mized effective potential (OEP) method 37[ or the local- 
ized Hartree-Fock (LHF) [30[ approximation. A recent 
work showed in fact that accurate photoelectron spec- 
tra can be calculated for very small gold clusters (up to 
4 atoms) by employing the LHF method together with 
a Lee-Yang-Paar (LYP) [38| correlation functional |39| . 
This approach is SIE- free and provides a good agreement 
with reference experimental data without the need for a 
rescaling of the orbital energies. It is thus potentially a 
valuable tool for the simulation of photoelectron spectra 
and the interpretation of PES experiments. However, de- 
spite the LHFLYP approach is much more accurate than 
conventional GGA methods, it has still two drawbacks: 
i) the correlation part of the functional is treated at the 
GGA level, and ii) it is yet a KS method, thus based on 
single-particle orbitals: even with the exact (unknown) 
XC functional the KS energy-level does not correspond 
to the exact PES levels [40[ , the only exception being the 
highest occupied molecular orbital (HOMO), which can 
be shown to be exactly related to the IP |4CH42l | . 

Alternative DFT approaches for the description of the 
PES levels are based on hybrid functionals, where a frac- 
tion of the non-local Hartree-Fock exchange is included 
in the Hamiltonian. This approach goes beyond the KS 
scheme and it is commonly applied for organic molecules. 
However, it has several problems when applied to metal- 



lie systems, because of the drawbacks of the Hartree- 
Fock method in this case |43l445| . Therefore its use is 
essentiaUy Umited to small clusters. Moreover, hybrid 
approaches are mostly based on an error cancellation ef- 
fects and thus can neither fully correct self-interaction er- 
ror nor can provide systematic improvements, while they 
introduce several additional parameters. Reliable results 
can be obtained instead using range separated hybrid 
functionals [46[, which take advantage of a screening of 
the Hartree-Fock exchange by separating short- and long- 
range effects. 

The fundamental solution to probe PES experiments 
dire ctly is to compute the quasiparticle energy spectrum 
471 ksf , as described by the quasiparticle equation 



i7o(r)0r (r)+/ dr'l].e(r, r', 6^ )<^r(r') = 



QPaQP 



(1) 

where iJo = — (1/2)V^ -I- Vext + Vh is the single-particle 
Hamiltonian, Vext is the external (nuclear) potential, Vh 
is the Hartree potential and Yjxc is the self-energy term: 
a non-local, non Hcrmitian and energy dependent oper- 
ator resulting from all electron many-body interactions. 
Here e^ and ip^ are the quasiparticle (QP) energies 
and eigcnstatcs, respectively. Unfortunately, calculations 
of the exact self energy for real systems are computa- 
tionally infeasible and approximations are required. We 
have calculated the self-energy operator within the GW 
approximation [43, [50| . 



^r,r',.) = - 



e'"'"^G'(r, 



,uj-uj')W{v,Y',uj')duj' 



(2) 

where G is the single particle Green's function and W is 
the dynamically screened Coulomb potential. 

In principle, Eq. ^ requires a self-consistent solution 
63|,|64|. Nevertheless, for many systems the DFT wave 
functions arc good approximations for the QP wavefunc- 
tions [30J, i40, j48j and the quasiparticle energies ej^^ can 
be obtained pertubatively ("one-shot" or GqWo calcula- 
tions) from the starting one-electron energies e" as: 



.QP 
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Originally, Go Wo has been applied in solids and, recently, 
calculations on bulk gold showed a QP band structure in 
good agreement with experimental PES [5l|- Moreover, 
this approach has been used also for accurate calculations 
of the IP of finite-size systems, both in closed-shell 52- 
73j and open-shell [52|, [7J| species. 

In this paper, we aim at assessing the merits and lim- 
itations of the different methods illustrated above and 
perform a comparative study of various computational 
approaches for the simulation of the first PES peak (i.e. 
IP) of several anionic gold clusters (up to 11 atoms; see 
fig. [1} for which experimental data are available. 






Au„(a) 




Aun(b) 



FIG. 1: Structures of gold cluster anions 



COMPUTATIONAL DETAILS 



Reference IP energies and the general structure of the 
clusters were taken from Ref. [12[. Geometries were 



optimized using the PBEint XC functional 75|, which 



has been shown to yield accurate results for such sys- 
tems |76|. For the larger clusters (Auf , Au^ and Au7i) 
we considered the lowest two structures from Ref. jl2| . 
However, for Auf at the PBEint/def2-TZVP level both 
isomers converged to the same geometry (Auf A in Ref. 
13). This was also confirmed by TPSS/dcf2-QZVPP 
calculations 771, l78[ , which failed to identify a minimum- 
energy structure corresponding to Auf B. Thus, only one 
structure for Au^ is considered here. The IP of differ- 
ent systems was obtained from the HOMO energy of the 
following methods: Hartree-Fock (HF), LHF 0, LH- 
FLYP [ifllii, PBEint [lil, PBEO ^. In addition, we 
considered the A SCF approach (PBEint functional) and 
the post-DFT many-body perturbation theory GqWo ap- 
proach (from PBEint orbitals; GoWo@PBEint). Finally, 
to asses the role of the reference orbitals for GqWo cal- 
culations we performed our perturbative treatment also 
starting from self-consitent static GW one-electron en- 
ergy and orbital results (GqWoOCOHSEX), which have 
been recently shown to yield in solid-state physics sim- 
ilar results as the more computationally expensive self- 
consistent GW calculations 80|, |81| and good agrement 
with the experimental transport properties of strongly 



correlated systems [6J, 1811 ISl The COHSEX [801 is in 



fact a static version of the GW approximation, where 
the self-energy can be written as 



I](r,r') = -^0?^(r)0?^(r')W^(r,r',c. = O) 



+ -^(r - y'){W{v, r', w = 0) - w(r - r')) (4) 



i.e. as a sum of an instantaneous screened exchange 
(SEX) term and a Coulomb hole (COH) correlation term. 
All DFT calculations were performed with the TUR- 
BOMOLEprogram package $^, using the def2-TZVP 
basis set [78[ and a scalar relativistic effective core po- 
tential (ECP) replacing the inner 60 electrons of each 
gold atom [8^. The GqWo and the self-consistent COH- 
SEX calculations were performed by using the FHI-AIMS 
program [8^, |8g| with the fourth-tier tight basis set (in- 
cluding additional /- and g-type functions) and scalar 
relativistic effects included through the zeroth-order rel- 
ativistic approximation (ZORA) [87|. With this ap- 
proach we estimate the GqWo orbital energies to be con- 
verged within 0.05 eV. Test calculations for the PBEO 
and PBEint functionals showed that the results obtained 
by the two programs also differ by less than 0.05 eV. In 
GW calculations the self energy E was first calculated 
on the immaginary frequency w-grid and then analyti- 
cally continued to the real axis grid as implemented in 
FHI-AIMS code 



RESULTS AND DISCUSSION 

The calculated IPs, as obtained with different meth- 
ods, are listed in Tab. HI together with the experimen- 
tal results. For Aug" and Au^j^ we report two isomers, 
because they are rather close in energy (at the PBEint 
level the (b) isomers are 0.21 and 0.07 eV higher for Aug" 
and Auj";^, respectively). For Aug (b) and Au^;^ COH- 
SEX results are not present, because of the excessive 
computational cost. The last line of Tab. U reports the 
mean absolute error (MAE) with respect to experiment. 
This is calculated including all clusters but Au^(b) and 
Auj~j^(a), because for these two clusters our best estimates 
(ASCF and GoWo@PBEint) show that the IP is signifi- 
cantly lower than the experimental one, suggesting that 
these clusters are not measured in experiments. 

The results confirm that the HOMO energy from GGA 
DFT calculations (PBEint results) is a very poor estima- 
tion of the IP, because the SIE leads to a strong under- 
estimation of the orbital energies. Much better results 
are found instead when the Kohn-Sham calculations are 
performed using an effective exact-exchange approach, as 
the LHF one. In this case in fact the MAE drops to 0.66 
eV and even further improvement can be achieved by 
inclusion of explicit correlation contributions (LHFLYP 
has a MAE of only 0.45 eV in line with the findings of 
Ref. [39[). We note, nonetheless, that this improvement 
is not systematic, as shown by the fact that LHFLYP 
performs well for the smaller clusters but it is worst than 
LHF for the largest ones. A possible reason for such a be- 
havior traces back to the limitations of the semilocal LYP 
correlation potential [32| (similarly to other semilocal po- 
tentials). It is not fully compatible with exact-exchange 
or OEP approaches [89|, and, more important, fails to 



properly describe the correlation effects in extended sys- 
tems of metallic character [90| . The net effect of the LYP 
functional is always to increase (in absolute value) the 
HOMO energy by approximately 0.7 eV, independently 
on the system under consideration. On the other hand, 
GqWo calculations suggest that the correlation effects are 
twice as big for the Au~ atom as for the largest clusters 
considered in this work, in agreement with the LHF vs. 
experimental trend. 

The importance of correlation effects is also remarked 
by a comparison of the LHF results with the Hartree- 
Fock ones, as the latter display a larger overall MAE. 
The difference can be in fact attributed not only to the 
local KS approximation (fuU-OEP exact-exchange Kohn- 
Sham results are often closer to HF than LHF ones [30[), 
but also to the LHF approximations to the full OEP, 
which make effective correlation contributions appear in 
the KS potential, so that an effective partial screening of 
the exchange is obtained [39[. On the contrary, in the 
Hartree-Fock method the correlation is completely lack- 
ing (absence of any screening) . Thus, Hartree-Fock calcu- 
lations yield an overestimation of the electron-electron re- 
pulsion in metallic systems and the IP obtained through 
the HOMO energy is underestimated (for the same reason 
the s-d splitting is overestimated in Hartree-Fock calcu- 
lations on metallic systems 9l|). This effect is instead 
mitigated in LHF calculations which in fact yield larger 
values for the IP. Finally, poor results are also obtained 
by considering the HOMO energy from hybrid calcula- 
tions, as both GGA and HF results underestimate the 
reference experimental values, so that no benefit from 
error compensation can be obtained. 

It is thus clear that to achieve a HOMO energy of 
high accuracy a proper treatment of exchange and cor- 
relation is necessary. This is indeed obtained within the 
GoWo(@PBEint) approach which in fact yields a MAE of 
only 0.14 eV, similarly with the ASCF approach which 
is based on total energy difi^erences (even slightly bet- 
ter Go Wo results are found using PBE starting orbitals, 
with a MAE of 0.11 eV; note however that the difference 
is below the estimated accuracy for Go Wo calculations). 
We remark that the ASCF method can be only applied 
for the HOMO energy and few other orbital energies cor- 
responding to the highest of each irreducible symmetry 
representation. On the contrary, the GqWo can be ap- 
plied straightforwardly for the calculation of all the QP 
energies. We note however that the perturbative Go Wo 
success in these calculations is not trivial: this method, 
originally applied to periodic systems, has quite few stud- 
ies about molecules and clusters where the inhomogenity 
of the density and low coordination number of atoms play 
an important role. 

In this sense it is interesting to observe that for the 
present gold cluster Go Wo calculations on top of self- 
consistent COHSEX (GW@COHSEX) show rather worse 
results than the GW@PBEint ones, in contrast to what 



TABLE I: Ionization potential for several anionic gold clusters as resulting from different a ppr oaches: HF, LHF, LH- 
FLYP,PBEint, ASCF, PBEO, COHSEX, GoWo@COHSEX, GoWo@PBEint. Experimental data [g are also reported. The 
last line shows the mean absolute error (MAE), excluding Au^(b) and Au^i(a). All results are in eV. 





HF 


LHF LHFLYP PBEint ASCF PBEO COHSEX " 


GoWo 






Cluster 


©COHSEX @ 


PBEint 


Exp. 


Au- 


0.80 


0.85 


1.51 


-0.84 


2.09 


-0.19 


3.00 


1.91 


2.23 


2.33 


Aug- 


2.46 


2.88 


3.57 


1.20 


3.57 


1.82 


4.84 


3.69 


3.92 


3.88 


Au^ 


1.91 


2.33 


3.02 


1.01 


3.03 


1.43 


4.08 


2.84 


2.87 


3.09 


Auf 


2.15 


2.43 


3.14 


1.59 


3.39 


1.89 


4.45 


3.25 


3.30 


3.46 


Au^(a) 


2.89 


3.55 


4.28 


1.95 


3.71 


2.45 


4.98 


3.75 


3.57 


3.83 


Aug-(b) 


2.36 


2.83 


3.56 


1.73 


3.44 


2.09 


- 


- 


3.30 


(3.83) 


Au-(a) 


2.69 


3.23 


3.96 


2.03 


3.65 


2.39 


- 


- 


3.52 


(3.80) 


Au-,(b) 


3.05 


3.77 


4.50 


2.18 


3.81 


2.62 


- 


- 


3.73 


3.80 


MAE 


1.19 


0.76 


0.45 


2.22 


0.14 


1.73 


0.95 


0.23 


0.14 





observed in solid-state calculations [81|. This is because 
of the poor quality of the COHSEX orbitals which over- 
estimate quite significantly the reference ones. This be- 
havior of the orbitals may trace back to the fact that in 
all clusters the HOMO is essentially described by a linear 
combination of s and p atomic orbitals, with negligible 
contribution from d and / atomic orbitals. In fact, scf- 
COHSEX was demonstrated to be a good starting point 
for Go Wo when localized d and / orbitals 8l| are im- 
portant (e.g., in strongly correlated systems). However, 
it can be expected to yield a poor description of delocal- 
ized electrons as those lying in s and p orbitals, due to 
the lack of dynamical screening. 



Finally, it is worthwhile to note that for the present 
case of gold anionic clusters the Hartree-Fock calcula- 
tions always underestimate the reference IP, whereas the 
Go Wo (and COHSEX) calculations yield an increase (in 
absolute value) of the HOMO energy with respect to HF. 
This behavior is opposite with respect to that usually ob- 
served in molecular systems where HF overestimates the 
IP due to the too negative value of the non-local exchange 
and GW has the effect of reducing the HOMO energy 
thanks to the screening of the exchange 1631 1. However, 
many exceptions exist, e. 
[5J1 and metal clusters 



in atoms [5^ |70[ , molecules 
62l . I63l l73{ . This latter behavior 



may be related to the important role of the correlation 
in systems such as the gold clusters studied here, which 
overtakes the screening. In fact, COHSEX calculations 
for the Au^ atom show that the effect of the screening 
term is to shift-up the HOMO energy by 0.2 eV, whereas 
the Coulomb hole shifts it down by 2.1 eV. The anionic 
character of the clusters instead appears to play a neg- 
ligible role, as similar findings are found also for neutral 
gold clusters |73j . 



CONCLUSIONS 

In conclusion, we have shown that refined many body 
calculations at the Go Wo ©PBEint level can describe the 
ionization potential of gold clusters with an accuracy 
of about 0.1 eV, outperforming any other orbital-based 
method. This result is a non-trivial confirmation of the 
goodness of the GW approach for metallic clusters and 
shows that the attainable accuracy is sufficient to dis- 
tinguish between different isomeric forms of a cluster of 
a given size. Thus the GoWq approach can be used to 
select specific clusters geometries by comparison with ex- 
perimental measurements (as in the case of Aug and 
Au]"]^). Similar results can be achieved at the ASCF level, 
but losing any orbital picture. Finally, within the Kohn- 
Sham scheme, relatively accurate results can also be ob- 
tained when the effective exact exchange LHF method 
is employed (eventually in conjunction with a semilocal 
correlation functional). This approach has a consider- 
ably lower computational cost than the GqWo and thus 
may be considered as a cheaper alternative to the latter. 



ACKNOWLEDGMENTS 

This work was partially funded by the European Re- 
search Council (ERG) Starting Grant FP7 Project DE- 
DOM, grant agreement no. 207441. We thank TURBO- 
MOLE GmbH for providing the TURBOMOLE program 
package and M. Margarito for technical support. 



[1] H. Schmidbaur, Gold Bulletin 23, 11 (1990). 

[2] M.-C. Daniel, D. Astruc, Chem. Rev. 104, 293 (2004). 

[3] P. Pyykko, Angew. Chem. Int. Ed. 43, 4412 (2004). 



[4] P. Pyykko, Inorganica Chimica Acta 358, 4113 (2005). [36 

[5] P. Pyykko, Chem. Soc. Rev. 37, 1967 (2008). 

[6] W.-W. Yam, E. C.-C. Cheng, Chem. Soc. Rev. 37, 1806 [37 

(2008). 
[7] G. J. Hutchmgs, M. Brust, H. Schmidbaur, Chem. Soc. 

Rev. 37, 1759 (2008). 
[8] Gold Chemistry. Applications and Future Directions in [38 

the Life Sciences, edited by F. Mohr (Wiley- VCH, Wein- 

heim, 2009). [39 

[9] B. Hammer, J. K. Norskov, Nature 376, 238 (1995). 

[10] D. Schoos, P. Weis, O. Hampe, M. M. Kappes, Phil. [40 

Trans. R. Soc. A 368, 1915 (2010). 
[11] K. J. Taylor, C. L. Pettiettehall, O. Cheshnovsky, R. E. [41 

Smalley, J. Chem. Phys. 96, 3319 (1992). 
[12] H. Haakkinen, B. Yoon, U. Landman, X. Li, H. J. Zhai, [42 

L. S. Wang, J. Phys. Chem. A 107, 6168 (2003). [43 

[13] J. Li, X. Li, H. J. Zhai, L. S. Wang, Science 299, 864 

(2003). [44 

[14] S. Bulusu, X. Li, L. S. Wang, X. C. Zeng, Proc. Natl. 

Acad. Sci. U.S.A. 103, 8326 (2006). [45 

[15] B. Yoon, P. Koskinen, B. Huber, O. Kostko, B. 

von Issendorff, H. Hkkinen, M. Moseler, U. Landman, [46 

ChemPhysChem 8, 157 (2007). 
[16] H. Hakkinen, Chem. Rev. 37, 1847 (2008). [47 

[17] N. Shao, W. Huang, Y. Gao, L.-M. Wang, X. Li, L. S. 

Wang, X. C. Zeng, J. Am. Chem. Soc. 132, 6596 (2010). [ 

[18] W. Huang, R. Pal, L.-M. Wang, X. C. Zeng, L. S. Wang, 

J. Chem. Phys. 132, 054305 (2010). [49 

[19] F. Furche, R. Ahlrichs, P. Weis, C. Jacob, S. Gilb, T. [50 

Bierweiler, M. M. Kappes, J. Chem. Phys. 117, 6982 

(2002). [51 

[20] P. Gruene, D. M. Rayner, B. Redlich, A. F. G. van der 

Meer, J. T. Lyon, G. Meijer, A. Fielicke, Science 321, 

674 (2008). [52 

[21] A. Fielicke, A. Kirilyuk, C. Ratsch, J. Behler, M. Schef- 

fler, G. von Helden, G. Meijer, Phys. Rev. Lett. 93, [53 

023401 (2004). [54 

[22] X. Xing, B. Yoon, U. Landman, J. H. Parks, Phys. Rev. 

B 74, 165423 (2006). [55 

[23] A. Lechtken, C. Neiss, J. Stairs, D. Schoss, J. Chem. 

Phys. 129, 154304 (2008). [56 

[24] M. P. Johansson, A. Lechtken, D. Schooss, M. M. [57 

Kappes, F. Furche, Phys. Rev. A 77, 053202 (2008). 
[25] A. Lechtken, C. Neiss, M. M. Kappes, D. Schooss, Phys. [58 

Chem. Chem. Phys. 11, 4344 (2009). 
[26] A. Schweizer, J. M. Weber, S. Gilb, H. Schneider, D. [59 

Schooss, M. M. Kappes, J. Chem. Phys. 119, 3699 

(2003). [60 

[27] S. Gilb, K. Jacobsen, D. Schooss, F. Furche, R. Ahlrichs, 

M. M. Kappes, J. Chem. Phys. 121, 4619 (1994). [61 

[28] A. N. Gloess, H. Schneider, J. M. Weber, M. M. Kappes, 

J. Chem. Phys. 128, 114312 (2008). [62 

[29] W. Kohn, L. J. Sham, Phys. Rev. 140, A1133 (1965). [63 

[30] F. Delia Sala, A. Gorhng, J. Chem. Phys. 115, 5718 

(2001). [64 

[31] M. Lundberg, P. E. M. Siegbahn, J. Chem. Phys. 122, 

224103 (2005). [65 

[32] E. Fabiano, F. Delia Sala, J. Chem. Phys. 126, 214102 

(2007). [66 

[33] A. J. Cohen, P. Mori-Sanchez, W. Yang, Science 321, [67 

792 (2008). 
[34] F. Delia Sala, Theor. Chem. Ace. 117, 981 (2007). [68] 

[35] R. Stowasser, R. Hoifmann, J. Am. Chem. Soc. 121, 3414 

(1999). 



A. Nagy, H. Adachi, J. Phys. B: At. Mol. Opt. Phys. 33, 
L585 (2000). 

F. Delia Sala in Orbital- Dependent Exact-Exchange 
Methods in Density Functional Theory and Chemical 
Modelling Applications and Theory Vol. 7 (Royal Soci- 
ety of Chemistry, RSC Publishing, 2010), pp. 115-161. 

C. Lee, W. Yang, R. G. Parr, Phys. Rev. B 37, 785 
(1988). 

E. Fabiano, M. Piacenza, F. Delia Sala, Phys. Chem. 
Chem. Phys. 11, 9160 (2009). 

D. P. Chong, O. V. Gritsenko, E. J. Baerends, J. Chem. 
Phys. 116, 1760 (2002). 

J. P. Perdew, R. G. Parr, M. Levy, J. L. Balduz, Phys. 

Rev. Lett. 49, 1691 (1982). 

J. P. Perdew, M. Levy, Phys. Rev. B 56, 16021 (1997). 

E. Fabiano, L. A. Constantin, F. Delia Sala, Int. J. 
Quant. Chem. (2013), doi: 10.1002/qua.24042. 

R. P. Messmer, T. C. Caves, C. M. Kao, Chem. Phys. 
Lett. 90, 296 (1982). 

D. Post, E. J. Baerends, Chem. Phys. Lett. 86, 176 
(1982). 

T. Stein, H. Eisenberg, L. Kronik, R. Baer, Phys Rev. 
Lett. 105, 266802 (2010). 

G. Onida, L. Reining, A. Rubio, Rev. Mod. Phys. 74, 
601 (2002). 

P. Duffy, D. P. Chong, M. E. Casida, D. R. Salahub, 

Phys. Rev. A 50, 4707 (1994). 

L. Hedin, Phys. Rev. 139, A796 (1965). 

F. Aryasetiawan, O. Gunnarsson, Rep. Prog. Phys. 61, 
237 (1998). 

T. Rangel, D. Kecik, P. E. Trevisanutto, G. M. Rig- 
nanese, H. V. Swygenhoven, V. Olevano, Phys. Rev. B 
86, 125125 (2012). 

E. L. Shirley, R. M. Martin, Phys. Rev. B 47, 15404 
(1993). 

M. Rohlfing, Int. J. Quantum Chem. 80, 807 (2000). 

J. C. Grossman, M. Rohlfing, L. Mitas, S. G. Louie, M. 

L. Cohen, Phys. Rev. Lett. 86, 472 (2001). 

S. Ishii, K. Ohno, Y. Kawazoe, S. G. Louie, Phys. Rev. 

B 65, 245109 (2002). 

Y. Pavlyukh, W. Hiibner, Phys. Lett. A 327, 241 (2004). 

A. Stan, N. E. Dahlen, R. van Leeuwen, Europhys. Lett. 

76, 298 (2006). 

M. L. Tiago, J. R. Chelikowsky, Phys. Rev. B 73, 205334 

(2006). 

E. Kikuchi, S. Ishh, K. Ohno, Phys. Rev. B 74, 195410 
(2006). 

Y. Noguchi, S. Ishu, K. Ohno, T. Sasaki, J. Chem. Phys. 
129, 104104 (2008). 

M. L. Tiago, J. C. Idrobo, S. Ogiit, J. Jellinek, J. R. 
Chelikowsky, Phys. Rev. B 79, 155419 (2009). 

F. Bruneval, Phys. Rev. Lett. 103 176403, (2009). 

C. Rostgaard, K. W. Jacobsen, K. S. Thygesen, Phys. 

Rev. B 81, 085103 (2010). 

M. Strange, C. Rostgaard, H. Hakkinen, K. S. Thygesen, 

Phys. Rev. B 83, 115108 (2011). 

L. Chiodo, M. Salazar, A. H. Romero, S. Laricchia, F. 

Delia Sala, A. Rubio, J. Chem. Phys. 135, 244704 (2011). 

S. Ke, Phys. Rev. B 84, 205415 (2011). 

X. Blase, C. Attaccalite, V. Olevano, Phys. Rev. B 83, 

115103 (2011). 

C. Faber, C. Attaccalite, V. Olevano, E. Runge, X. Blase, 

Phys. Rev. B 83, 115123 (2011). 



[69] S. Sharifzadeh, I. Tamblyn, P. Doak, P. T. Darancet, J. 

B. Neaton, Eur. Phys. J. B 85, 323 (2012). 
[70] F. Bruneval, J. Chem Phys. 136, 194107 (2012). 
[71] P. Umari, S. Fabris, J. Chem Phys. 136, 174310 (2012). 
[72] N. Marom, F. Caruso, X. Ren, O. T. Hofmann, T. 

Korzd6rfer,J. R. Chelikowsky, A. Rubio, M. Scheffler, P. 

Rinke, Phys. Rev. B 86, 245227 (2012). 
[73] M. J. van Setten, F. Weigend, F. Evers, J. Chem. Theor. 

Comput. 9, 232 (2013). 
[74] J. Lischner, J. Deshppe, M. Jain, S. G. Louie, Phys. Rev. 

Lett. 109, 036406 (2012). 
[75] E. Fabiano, L. A. Constantin, F. Delia Sala, Phys. Rev. 

B 82, 113104 (2010). 
[76] E. Fabiano, L. A. Constantin, F. Delia Sala, J. Chem. 

Phys. 134, 194112 (2011). 
[77] J. Tao, J. P. Perdew, V. N. Staroverov, G. E. Scuseria, 

Phys Rev. Lett. 91, 146401 (2003). 
[78] F. Weigend, F. Furche, R. Ahlrichs, J. Chem. Phys. 119, 

12753 (2003). 
[79] C. Adamo, V. Barone, J. Chem. Phys. 110, 6158 (1999). 
[80] F. Bruneval, N. Vast, L. Reining, Phys. Rev. B 74, 

045102 (2006). 
[81] M. Gatti, F. Bruneval, V. Olevano, L. Reining, Phys. 

Rev. Lett. 99, 266402 (2007). 
[82] T. Rangel, A. Ferretti, P. E. Trevisanutto, V. Olevano, 



G.-M. Rignanese, Phys. Rev. B 84, 045426 (2011). 
[83] TURBOMOLE V6.4 2012, a development of University 

of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 

1989 -2007, TURBOMOLE GmbH , since 2007; available 

from |http : //www . turbomole . com[ 
[84] D. Andrae, U. Haussermann, M. Dolg, H. StoU, H. 

Preuss, Theor. Chim. Acta 77, 123 (1990). 
[85] V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. 

Ren, K. Renter, M. Scheffler, Comput. Phys. Commun. 

180, 2175 (2009). 
[86] V. Havu, V. Blum, P. Havu, M. Scheffler, J. Comput. 

Phys. 228, 8367 (2009). 
[87] S. Faas, J. G. Snijders, J. H. van Lenthe, E. van Lenthe, 

E. J. Baerends, Chem. Phys. Lett. 246, 632 (1995). 
[88] X. Ren, P. Rinke, V. Blum, J. Wieferink, A. Tkatchenko, 

A. Sanfilippo, K. Renter, M. Scheffler, New Journal of 

Physics 14, 053020 (2012). 
[89] L. A. Constantin, E. Fabiano, F. Delia Sala, Phys. Rev. 

B 86, 035130 (2012). 
[90] J. Paier, M. Marsman, G. Kresse, J. Chem. Phys. 127, 

024103 (2007). 
[91] F. Delia Sala, E. Fabiano, S. Laricchia, S. D'Agostino, 

M. Piacenza, Int. J. Quant. Chem. 110, 2162 (2010). 



